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It has been known for a long time that many simple liquids have surprisingly similar structure 
as quantified, e.g., by the radial distribution function. A much more recent realization is that the 
dynamics are also very similar for a number of systems with quite different pair potentials. Systems 
with such non-trivial similarities are generally referred to as “quasi-universal”. From the fact that 
the exponentially repulsive pair potential has strong virial potential-energy correlations in the low- 
temperature part of its thermodynamic phase diagram, we here show that a liquid is quasi-universal 
if its pair potential can be written approximately as a sum of exponential terms with numerically 
large prefactors. Based on evidence from the literature we moreover conjecture the converse, i.e., 
that quasi-universality only applies for systems with this property. 


INTRODUCTION 

A cornerstone of liquid-state theory is the fact that the hard-sphere (HS) model gives an excellent representation 
of simple liquids. This was discussed in 1959 by Bernal [TJ, and the first liquid-state computer simulations at the 
same time by Alder and Wainwright likewise studied the HS system |2J. Following this the HS model has been the 
fundament of liquid-state theory, in particular since the 1970s when perturbation theory matured into its present form 
[5H5]. The HS model is usually invoked to explain the fact that many simple liquids have very similar structure. This 
cannot explain, however, the more recent observation of dynamic quasi-universality (SHE! or why some mathematically 
simple pair potentials violate quasi-universality [151 - 119] . 

An early hint of dynamic quasi-universality was provided by Rosenfeld who in 1977 showed that the diffusion 
constant is an almost universal function of the excess entropy (the entropy minus that of an ideal gas at the same 
temperature and density) 0. This finding did not attract a great deal of attention at the time, but the last decade has 
seen renewed interest in excess-entropy scaling III HI and, more generally, in the striking similarities of the structure 
and dynamics of many simple model liquids mm- Thus Heyes and Branka and others have documented that inverse 
power-law (IPL) systems with different exponents have similar structure and dynamics [9i 1TT1120H251 . Medina-Noyola 
and coworkers have established that quasi-universality extends to the dynamics for Newtonian as well as Brownian 
equations of motion nfflnmna, Scopigno and co-workers |2l] have documented HS-like dynamics in liquid gallium 
studied by quasi-elastic neutron scattering, and Liu and coworkers have suggested a mapping of a soft-sphere system’s 
dynamics to that of the HS system [123 ■ 

Further examples of similarities between apparently quite different systems, identified throughout the years but 
still largely unexplained, include: the Young-Andersen approximate scaling principle according to which two systems 
that have the same radial distribution function - even at different thermodynamic state points - also have the same 
dynamics [24j : the fact that different systems have similar order-parameter maps in the sense of Debenedetti and 
coworkers 0125], i.e., that when the relevant orientational order parameter is plotted against the translational one, 
almost identical curves result; the Lindemann and other melting rules |26l l27j ; freezing rules like Andrade’s finding 
that freezing initiates when the reduced viscosity upon cooling reaches a certain value [28] or the Hansen-Verlet rule 
that a liquid freezes when the maximum structure factor reaches 2.85 [251 IBP] , 

In the literature the term “quasi-universality” refers to the general observation that different simple systems have 
very similar physics as regards structure and dynamics, but before proceeding we need to define the term precisely. To 
do this, recall first the definition of so-called reduced quantities. These are quantities made dimensionless by dividing 
by the appropriate combination of the following three units: the length unit is p~ 1,/3 where p is the number density, 
the energy unit is ksT where T is the temperature, and the time unit is p~ x ^ ^/m/ksT for Newtonian dynamics 
where m is the particle mass (for Brownian dynamics a different time unit is used [5ll 132] ). Henceforth, reduced 
quantities are denoted by a tilde, for instance the reduced pair distance is defined by f = p l ^r. 

By quasi-universality we shall mean the property of many simple model liquids that knowledge of a single quantity 
characterizing the structure or dynamics in reduced units is enough to determine all other reduced-unit structural and 
dynamic quantities to a good approximation. Note that a special case of this is Rosenfeld’s observation that the excess 
entropy over ks, a structural quantity, determines the reduced diffusion constant. Likewise, the Young-Andersen [55] 
approximate scaling principle is a consequence of quasi-universality as defined here - as these authors expressed it 
“... certain dynamical properties are very insensitive to large changes in the interatomic potential that leave the pair 
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correlation function largely unchanged”. 

It is straightforward to show that if quasi-universality applies, the other points mentioned above also follow. The 
fundamental questions are: Which systems are quasi-universal? What causes quasi-universality? The conventional 
explanation of quasi-universality starts from the fact that the HS system provides an excellent reference system for 
structure calculations ei m mm m- Approximate theories of liquid dynamic properties like renormalized kinetic 
theory or mode-coupling theory in its simplest version predict that the dynamics is uniquely determined by the static 
structure factor, i.e., that structure determines dynamics. This reasoning led to the search for and discovery of 
dynamic quasi-universality back in 2003 (refs. [TOl f24j> . 

The dynamics of the HS system consists of constant-velocity free-particle motion interrupted by infinitely fast 
collisions, which is quite different from the continuous motion described by Newton’s laws for smooth potentials. 
Thus it is far from physically obvious why the HS explanation of quasi-universal structure extends to the dynamics. 
Moreover, while the HS explanation does account for the finding that, e.g., the Gaussian-core model da nn ns 
violates quasi-universality - a model without harsh repulsions except at low temperatures - it cannot explain why 
some strongly repulsive or hard-core pair-potential systems violate quasi-universality. Examples of such systems are 
the Lennard-Jones Gaussian model m and the Jagla model m- Finally, the one-component plasma model is well 
known to be quasi-universal [6] , but it is not intuitively clear when and why the gently varying Coulomb force can be 
well approximated by the harsh HS interaction. 

In this paper we do not use the HS reference system. In a continuation of recent works [381 - 140] we take a different 
approach to quasi-universality by showing that any pair potential, which can be approximated by a sum of exponential 
terms with numerically large prefactors, is quasi-universal. This introduces the “EXP quasi-universality class” of pair 
potentials. Based on the available evidence from the literature we moreover conjecture the converse, that is, that all 
quasi-universal systems are in the EXP quasi-universality class. 


RESULTS 

The exponentially repulsive pair potential 

We study below the monatomic system described by the purely repulsive “EXP” pair potential 


Wxp (r;e,a) = ee r, ° . (1) 

This was discussed already in 1932 by Born and Meyer and by Buckingham in 1938 in the context of a pair potential 
with an exponentially repulsive term plus an r -6 attractive term EU 02]; note also that the well-known Morse pair 
potential is a difference of two exponentials [43]. The EXP pair potential may be justified physically as reflecting the 
overlap of electron wavefunctions in conjunction with the fact that bound-state wavefunctions decay exponentially 
in space jJTj. Although the EXP pair potential since 1932 has been used occasionally in computer simulations and 
for interpretation of experiments mmm, it never became a standard pair potential like the Lennard-Jones and 
Yukawa pair potentials [251 02 ] ■ Given the mathematical simplicity of the EXP function this may seem surprising, 
but a likely explanation is that no systems in nature are believed to be well described by this purely repulsive pair 
potential. The present paper suggests that it may nevertheless be a good idea to regard the EXP pair potential as 
the fundamental building block of simple liquids’ physics, much like the exponential function is in the Fourier and 
Laplace transform theories of pure mathematics. 

Our simulations focused on the low-temperature part of the thermodynamic phase diagram, i.e., where 

k B T « e . (2) 

The potential energy as a function of all the particle coordinates R = (iq, ...,rjv) is denoted by C/exp(R)- Figure 
1 shows scatter plots of [/exp(R) for several configurations R plotted versus the potential energies of the same 
configurations scaled uniformly to a different density, i.e., C/exp(R) plotted versus Lexp(AR) in which A 3 is the 
ratio of the two densities in question. The figure was constructed by simulating five temperatures at the density 
pc r 3 = 1.0 • 10 -3 (lower line, different colors); at each of these five state points configurations were selected and scaled 
to four higher densities. We see that there are very strong correlations between the potential energies of scaled and 
unsealed configurations. 
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FIG. 1: Scatter plots of potential energies per particle of the EXP system. At the state point with density pa 4 = 1.0- 10~‘ ! and 
temperature ksT/e = 1.0 • 1CP 4 several configurations were selected from an equilibrium simulation and scaled uniformly to 
four higher densities (the light brown points); the same was done for four other temperatures at this density. The figure shows 
that there are very strong correlations between scaled and unsealed potential energies, with a scaling factor that only depends 
on the two densities involved. The lines are best fits to the green data points, that is, those generated from simulations at 
temperature IcbT/s = 1.4 • 10~ 4 . This figure validates the hidden-scale-invariance property of the EXP pair-potential system 
as expressed below in Eq. Q. 


As shown elsewhere strong correlations between scaled and unsealed potential energies are characteristic for systems 
that have strong correlations between their virial and potential-energy thermal equilibrium fluctuations at the relevant 
state points m- Systems with this property include not only “atomic” pair-potential systems like the IPL, Yukawa, 
and Lennard-Jones (LJ) systems, etc, but also a number of rigid molecular model systems and even the flexible 
LJ chain model [501 [31]. Such systems were previously referred to as “strongly correlating”, but are now called 
Roskilde-simple or just Roskilde systems |52H57j which avoids confusion with strongly correlated quantum systems. 

The strength of the virial potential-energy correlations of the EXP system is reported in Fig. 2(a) for a large 
number of state points. If virial and potential energy are denoted by W and C/, respectively, and sharp brackets denote 
canonical constant-volume ( NVT) averages, the color coding of the figure gives the Pearson correlation coefficient R 

defined [55] by 


(AWAU) 

V<(acto(awtT 


(3) 


A pragmatic definition is that a given system is Roskilde-simple at the state point in question if R > 0.9 [55]. This is 
the case for the EXP system whenever ksT/e < 0.1. 

Figure 2(b) zooms in on the low-temperature part of the phase diagram where the EXP system has particularly 
strong virial potential-energy correlations. In both figures the line segment have slope 7 , the slope of the so-called 
isomorph through the state point in question that is calculated from the expression 7 = (AW AU) / ((AU) 2 ) [31] . 
An isomorph is a curve in the phase diagram along which structure and dynamics to a good approximation are 
invariant in reduced units [33[52]. A system has isomorphs if and only if the system has strong virial potential-energy 
correlations at the relevant state points m ■ The existence of isomorphs implies that the phase diagram becomes 
effectively one-dimensional for many physical properties; Ref. m gives a recent review of the isomorph theory with 
a focus on its validation in computer simulations and experiments. The curves in Fig. 2(a) are examples of isomorphs 
of which the full curve is the melting isomorph [31] . Isomorphs of the EXP system were studied briefly in Ref. 00] 

Before proceeding we address the question why the EXP pair potential has strong virial potential-energy correlations. 
Our explanation refers to IPL pair potentials, which by Euler’s theorem for homogeneous functions have 100% virial 
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FIG. 2: Thermodynamic phase diagram of the EXP system. Density and temperature are both given in the unit system defined 
by the pair-potential parameters (Eq. 0 ). (a) gives an overview of the phase diagram and (b) zooms in on its low-temperature 
part; in both cases the colors indicate the value of the virial potential-energy correlation coefficient R defined in Eq. The 
slope of each line segment is that of the isomorph through the state point in question, the curve along which structure and 
dynamics are (approximately) invariant in reduced units. The full curves mark the melting isomorph; this curve’s width in (a) 
is approximately that of the coexistence region. The dashed curves mark a few other isomorphs. 
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FIG. 3: Illustration of how the EXP pair potential (black curve) may be fitted by an eIPL pair potential [58]. The eIPL 
pair potential is a sum of an inverse power-law (IPL) term and a linear term (red dashed curve); f = rp 1 ' 3 is the reduced 
pair distance and the pair potential is given in units of ksT. The eIPL function was fitted to the EXP function over the 
range of radial distances for which the radial distribution function (RDF) at the state point defined by pa 3 = 1.0 • 10 -3 and 
ksT/e = 5.0 • 10~ 5 is larger than unity. The blue and green curves are the RDFs at this state point of the EXP pair potential 
and the IPL pair potential, respectively. The EXP and IPL pair potentials’ diffusion constants differ by less than 1% at this 
state point (data not shown). 


potential-energy correlations (recall that the microscopic virial is defined generally by W( R) = (—1/3)R • V£/(R) 
0)- Figure [ 3 ] refers to the state point pa 3 = 1.0 • 10 3 , ksT/e = 5.0 • 10 5 . The figure shows that the EXP pair 
potential (black curve) may be fitted very well over the entire first coordination shell by an extended IPL (“eIPL”) 
function, which is defined as an IPL term plus a linear term (the red dashed line). As shown in Ref. [58; the linear 
term contributes little to neither the virial nor the potential-energy constant-volume fluctuations because the sum 
of nearest-neighbor distances is almost constant. This is because if one particle is moved, some nearest-neighbor 
distances decrease and others increase, resulting in almost no change in the sum of the nearest-neighbor distances 
(this argument is exact in one dimension). As a consequence, the WU fluctuations are dominated by the IPL term 
and thus strongly correlating. This explains why the EXP pair potential has very strong WU correlations. To confirm 
the equivalence between the EXP pair potential and the IPL pair potential (1.13 f -8 ' 41 ), we performed a simulation 
with the IPL pair potential at the same state point. As predicted [59 EDI there is good agreement between the two 
systems’ structure (blue and green curves) and dynamics (not shown), confirming that the linear term of the eIPL 
pair potential fitting the EXP pair potential, 0.512 f — 0.738, is not important for the physics. 


Derivation of quasi-universality 

We proceed to show that the strong virial potential-energy correlation property of the EXP system implies quasi¬ 
universality for a large class of monatomic systems. As shown recently [39], any Roskilde-simple system, that is, 
system with strong virial potential-energy correlations and isomorphs, is characterized by “hidden scale invariance” 
in the sense that two functions of density exist, h(p) and g(p), such that the potential-energy function f7(R) can be 
expressed as follows 


C/(R)-/i(p)$(R)+ 5 (p). (4) 

Here $(R) is a dimensionless, state-point-independent function of the reduced dimensionless configuration vector 
R = y o 1 / 3 R. In particular, Eq. Q implies the strong correlations between the potential energies of scaled and unsealed 
configurations documented in Fig. 1: Consider two configurations Ri and R 2 at density p\ and p 2 , respectively, 
with the same reduced coordinates, i.e., obeying pJ^Ri = p^ 3 R 2 - By elimination of <f>(R) in Eq. Ml) we get 
£/(R 2 ) = {h{p 2 )/h(pi))U(Ri) + g(p 2 ) — g(pi)h(p 2 )/h(pi). Thus any two configurations with the same reduced 
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coordinates have potential energies that are approximately linear functions of each other with constants that only 
depend on the two densities in question, compare Fig. 1. 

According to Eq. Q a change of density implies that the potential-energy landscape to a good approximation 
undergoes a linear, affine transformation. Moreover, Eq. 0 implies invariant structure and dynamics in reduced 
units along the curves defined by constant h(p)/kBT, the equation identifying the isomorphs HUIMIEZI: along any 
curve defined by /i(p)/fcgT = C the reduced force F = —Vf7(R) /ksT is given by F = — C V$(R), which implies that 
F is a certain function of the reduced coordinates, i.e., F = F(R). Since Newton’s second law in reduced coordinates 
is d 2 R/dt 2 = F (the reduced mass is unity), it follows that the particles move in the same way at different isomorphic 
state points - except for the trivial linear uniform scalings of space and time involved in transforming back to real 
units. 

Since the EXP pair potential has strong virial potential-energy correlations, Eq. Q implies that functions h EXP (p), 
3bxp(p), and ^expCR) exist such that 


C^exp(R) — ^EXP (p)'f’EXp(R-) + 5exp(p) ■ (5) 

The functions h EXP (p) and g EXP (p) in Eq. ([5]) both have dimension energy. They can therefore be written h EXP (p) = 
£ h EXP (pcr 3 ) and g EXP (p) = eg EXP {pa 3 ) in which h EXP and g EXP are dimensionless functions that only depend on the 
dimensionless density pa 3 . 

Consider now a pair potential v(r) that can be expressed as follows 


v(r) = / f(a)e r ^da. (6) 

Jo 

The system’s potential energy is the sum of the individual pair-potential contributions with the same weights as 
in Eq. (JcJ) . Therefore, if U{ R) is the potential energy and we define h(p) = f(a)h EXP (pa 3 )da and g{p) = 
fo° f( cr )gEXp(p tj3 )da, it follows from Eq. (|HJ) that 


C/(R) = h(p)$EXp(R) + g{p) ■ (7) 

Because the reduced-unit physics is encoded in the function $exp(R) via Newton’s equation d 2 R/dt 2 = 
—C vI>exp(R) where C = h(p)/kBT identifes the isomorph through the state point in question, Eq. ([ 7 ]) implies 
identical structure and dynamics to a good approximation for systems obeying Eq. 0- This argument would be 
exact if the EXP pair potential had 100% virial potential-energy correlations. This is not the case, though, and the 
approximation Eq. ([7| is only useful when it primarily involves EXP functions from the low-temperature part of the 
phase diagram where strong virial potential-energy correlations and thus Eq. © apply. 

We have shown that a pair potential v{r) is quasi-universal if it can be written as a sum of low-temperature EXP 
pair potentials. To translate this into an operational criterion, note the following. If the integral in Eq. 0 is 
discretized into a finite sum and expressed in terms of the reduced pair potential v = v/ksT regarded as a function 
of the reduced pair distance (r = p 3 7 3 r), the condition for quasi-universality is that 


m = A,e-^ , |A,|»1. (8) 

i 

It is understood that the “wavevectors” Uj are not so closely spaced that large positive and negative neighboring terms 
may almost cancel one another. Several points should be noted: 1) Equation ^ is state-point dependent because the 
function v(f) varies with state point; 2) a continuous integral of EXP functions does not automatically obey Eq. ^ 
- it is necessary that the integral can be approximated by a finite sum of EXP terms, each with a numerically large 
prefactor; 3) a sum or product of two pair potentials obeying Eq. (|8]) with all A j > 0 gives a function that also obeys 
this equation. 


Important examples 

Consider first the IPL pair potential v n (r) = £(r/a) n . In terms of the reduced radius r, the reduced IPL pair 
potential is given by v n (r) = u n (r)/fcsT = T n f~ n in which T n = (per 3 )”/ 3 The mathematical identity 
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J 0 °° x n ~ x exp(— x)dx = (n — 1)! implies that v n (r) = [r„/(n — 1)!] J 0 °° u n_1 exp(— uf)du. Discretizing the integral 
leads to v n (f) = [r n /(n — 1)!] Au YlJLo (0 + l/2)A«) n_1 exp (—(J + 1/2)A uf). Writing ((j + l/2)Au) n_1 = exp[(n — 
1) ln((j + 1/2)Ait)], by differentiation with respect to j it is easy to see that the dominant contributions to the sum 
come from the terms with (n — 1 )/{j + 1/2) ~ A uf. Thus for typical nearest-neighbor distances (r ~ 1) the terms 
with (j + l/2)Au ~ n — 1 are the most important ones. For these values of j the prefactor of the exponential in 
the above sum is roughly T n Au(n — l) n-1 /(n — 1)!. The largest realistic discretization step A u is of order unity, 
so we conclude that for values of n larger than three or four Eq. (8|//s obeyed unless T n is very small, a condition 
that applies for the state points that have typically been studied 191 1201123] . The case of a Coulomb repulsive system 
(n = 1) is discussed in the next section. 

As a consequence of the above, at most state points the potential energy of the IPL system, C/ n (R), can be written 


Cn(R) = ^n(/3)^EXp(R) + 9n(p) ■ 


(9) 


Since t/ n (R) oc p n / 3 for the density variation induced by a uniform scaling of a configuration R, i.e., keeping R 
constant, one has h n (p) oc p n / 3 and g n (p) oc p n / 3 . This means that two numbers a n and f3 n exist such that 


C„(R) £* £ a„(pir 3 ) n / 3 l> E xp(R) + Pn(pv 3 ) n/3 


. For a general pair potential of the form v(r) = £ v n (r/a) n , by a 


linear combination of Eq. |9j) we arrive at Eq. ( 7| in which h(p) = £ J2 n v nOt n {p^ 3 ) n ^ 3 and g{p) = £ v nPn(po' 3 ) n ^ 3 ■ 

A well-known case is the Lennard-Jones (LJ) pair potential ulj(t) = 4£[(r/cr) -12 — (r/cr) -6 ]. Consider the LJ liquid 
state point given by per 3 = 1, fcgT = 2£. In terms of the above-defined reduced IPL functions, it is easy to see that 
this is of the form Eq. Q, implying quasi-universality of the LJ liquid at this state point. 

An application of the EXP theory of quasi-universality is the intriguing “additivity of melting temperatures” |6Tj 
according to which if two systems have melting temperatures that as functions of density are denoted by T mi i(p) and 
T m>2 (p), the melting temperature of the system with the sum potential energy is T m .i(p) + T m ^{p)- This property 
follows from quasi-universality because the dynamics of all three systems are controlled by the same function 'Fexp (R) • 
For the EXP system melting initiates when this function’s average upon heating reaches a certain value, and the same 
must apply for all quasi-universal systems. In particular, since for an IPL system one has T m (p) cx p ra//3 [62] , the 
melting temperature of the LJ system varies with density according to the expression T m (p) = Ap 4 — Bp 2 (refs. 

mm- 


DISCUSSION 

Denominating the systems that obey Eq. ([8]) collectively as the EXP quasi-universality class, we have shown that 
all systems in this class are quasi-universal in the sense of this paper: if a single reduced-unit structural or dynamic 
quantity is known, all other reduced-unit structural or dynamic quantities are known to a good approximation. This 
is because these are all encoded in the function $ E xp(R), and any reduced-unit quantity characterizing structure 
or dynamics identifies the constant C = h(p)/kBT of the reduced-unit version of Newton’s second law d 2 R/dt 2 = 
—C Vl> E xp(R). 

The obvious question is whether all quasi-universal systems are in the EXP quasi-universality class. We cannot 
prove this, but conjecture it is the case based on the available evidence in the literature for the following systems: 

The Jagla pair potential. This is a HS potential plus a finite-width potential well defined by two terms that are 
linear in r m • This pair potential is not in the EXP quasi-universality class - it cannot be approximated as a sum 
of exponentials because its Laplace transform only has a pole at zero. Indeed, the Jagla pair potential reproduces 
water’s anomalous density maximum and has a liquid-liquid critical point m, properties which are both inconsistent 
with quasi-universality. 

The Gaussian core model (GCM). This is a Gaussian centered at r = 0. For this system there is a reentrant 
BCC phase above the triple point m and the transport coefficients have a non-monotonous density dependence 
at constant temperature pT5] . These observations both contradict quasi-universality. Consistent with the above 
conjecture, a representation of the form exp(— r 2 /2a 2 ) = J 0 °° (j){u)e~ ur du does not exist because the Laplace transform 
of a Gaussian has no poles; thus the GCM system is not in the EXP quasi-universality class. At low temperatures 
and low densities the GCM pair potential can be approximated well by an exponential, however, and in this part of 
the phase diagram the GCM system is indeed quasi-universal [lTJ 15]. 

The LJ Gaussian (LJG) model. This pair potential is arrived at by adding a negative, displaced Gaussian to the LJ 
pair potential PH- Like the GCM it cannot be written as a sum of exponentials. The LJG system has thermodynamic 






and dynamic anomalies m and “a surprising variety of crystals” mi. both of which are observations that violate 
quasi-universality. 

The one-component plasma (OCP). This is the common name for the single-charge Coulomb system. The Coulomb 
potential is too long ranged for a thermodynamic limit to exist for the free energy per particle unless a uniform charge- 
compensating background is introduced nans]. Nevertheless, any finite OCP system is well defined and amenable to 
computer simulation in fact the OCP system was an important example in Rosenfeld’s original paper on excess- 
entropy scaling [5]. The OCP system is the n = 1 case of the above discussed IPL pair potential; it is quasi-universal 
in the dense fluid case, i.e., whenever Pi 1 (writing v\ = T\/r). Violations of quasi-universality are indeed known 
to gradually appear when Tj goes below 50 [53] , 

The Yukawa pair potential. This is given by the expression v(r) = e(u/r) exp(— r/a). Since v{r)/e = / 0 °° exp(— r(u+ 
1 )/a)du it is easy to see that this pair potential is in the EXP quasi-universality class whenever k^T -C e. This is 
consistent with known properties |65| . 

This paper has argued that one may replace the HS by the EXP system as the generic simple pair-potential system 
from which quasi-universality is derived. We do not wish to suggest that the HS system’s role in liquid-state physics is 
entirely undeserved, however, because this model is still uniquely simple and mathematically beautiful. Nevertheless, 
by using the EXP pair potential as the fundamental building block, a number of advantages are obtained. First of all, 
smoothness is ensured. Secondly, if the conjecture that all quasi-universal functions are in the EXP quasi-universality 
class is confirmed, we now have a mathematically precise characterization of all quasi-universal pair potentials (Eq. 
([8])). Thirdly, the present explanation of quasi-universality gives a natural explanation of dynamic quasi-universality. 
Finally, the role of the HS pair potential in liquid-state physics is clarified: quasi-universality is not caused by a given 
system’s similarity to the HS system - rather, quasi-universality implies similarity between the properties of many 
systems and those of the HS system simply because the latter is in the EXP quasi-universality class by being the 
n —» oo limit of an IPL system. - Note that while this paper focused on quasi-universality for liquid models, the 
crystalline phases of these models are likewise quasi-universal, for instance by having quasi-universal radial distribution 
functions, phonon spectra, and vacancy jump dynamics |66l . 

In continuing work it will be interesting to investigate the consequences of the EXP approach for the fact that 
quasi-universality appears to apply beyond the framework of thermal equilibrium, e.g., for the jamming transition 
in which the HS model is presently used as the generic model m ■ It is an open question whether replacing the HS 
model by the EXP model in thermodynamic perturbation theory would have advantages [51]. Another open question is 
to which extent mixtures are quasi-universal [68] . Finally, recent results from the works of Truskett and collaborators 
El show that classical Rosenfeld excess-entropy scaling may be modified into more general excess-entropy scalings, 
and it would be interesting to investigate whether the present approach can somehow be extended to account for this. 


METHODS 

A system of N = 1, 000 particles interacting via the EXP pair potential was simulated using standard Nose-Hoover 
NVT simulations with a time step of 0.0025 and thermostat relaxation time of 0.2 (LJ units). A shifted-forces 
cutoff at r = 2p~ x l^ was used for densities below 1.0 • 10 —3 , at larger densities the cutoff was 4p -1 / 3 . At each 
state point the simulations involved 10,000,000 time steps after equilibration. For temperatures below 2.0 • 10 -3 
the simulations were initiated from a state of 2,000 particles placed in a body-centered cubic crystal structure. 
The melting isomorph was determined by the interface pinning method |69| - the NPT simulations involved here 
were made using LAMMPS EH (http://lammps.s andia.gov) with shifted-potential cutoffs ranging from 2.5 p 1 ^ 3 to 
4.0 p -1 / 3 for a system of 2, 560 particles; the time step was 0.005 and the relaxation time was 0.4 for the thermostat 
and 0.8 for the barostat. 
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